The RHMC Algorithm for 2 Flavours of Dynamical Staggered Fermions

We describe an implementation of the Rational Hybrid Monte Carlo (RHMC) algorithm for dynamical computations with two flavours of staggered quarks. We discuss several variants of the method, the performance and possible sources of error for each of them, and we compare the performance and results to the inexact R algorithm.Comment: Lattice2003(machine) 3 pages, 1 figure. Added referenc

Testing Algorithms for Finite Temperature Lattice QCD

We discuss recent algorithmic improvements in simulating finite temperature QCD on a lattice. In particular, the Rational Hybrid Monte Carlo(RHMC) algorithm is employed to generate lattice configurations for 2+1 flavor QCD. Unlike the Hybrid R Algorithm, RHMC is reversible, admitting a Metropolis accept/reject step that eliminates the $\mathcal{O}(\delta t^2)$ errors inherent in the R Algorithm. We also employ several algorithmic speed-ups, including multiple time scales, the use of a more efficient numerical integrator, and Hasenbusch pre-conditioning of the fermion force.Comment: 4 pages, 2 figures, poster presented at International Conference on Strong and Electroweak Matter 2006 (SEWM2006), BNL, May 10-13, 200

Comparing the R algorithm and RHMC for staggered fermions

The R algorithm is widely used for simulating two flavours of dynamical staggered fermions. We give a simple proof that the algorithm converges to the desired probability distribution to within O(dt^2) errors, but show that the relevant expansion parameter is (dt/m)^2, m being the quark mass. The Rational Hybrid Monte Carlo (RHMC) algorithm provides an exact (i.e., has no step size errors) alternative for simulating the square root of the staggered Dirac operator. We propose using it to test the validity of the R algorithm for simulations carried out with dt m.Comment: 3 pages, proceedings from Lattice 2002 poster presentatio

The P Cygni supergiant [OMN2000] LS1 – implications for the star formation history of W51

Original article can be found at: http://www.aanda.org/ Copyright The European Southern Observatory (ESO) DOI: 10.1051/0004-6361/200911980Aims. We investigate the nature of the massive star [OMN2000] LS1 and use these results to constrain the history of star formation within the host complex W51. Methods. We utilised a combination of near-IR spectroscopy and non-LTE model atmosphere analysis to derive the physical properties of [OMN2000] LS1 , and a combination of theoretical evolutionary calculations and Monte Carlo simulations to apply limits on the star formation history of W51. Results. We find the spectrum of [OMN2000] LS1 to be consistent with that of a P Cygni supergiant. With a temperature in the range of 13.2–13.7 kK and log( ) , it is significantly cooler, less luminous, and less massive than proposed by previous authors. The presence of such a star within W51 shows that star formation has been underway for at least 3 Myr, while the formation of massive O stars is still on going. The lack of a population of evolved red supergiants within the complex shows that the rate of formation of young massive clusters at ages 9 Myr was lower than currently observed. We find no evidence of internally triggered, sequential star formation within W51, and favour the suggestion that star formation has proceeded at multiple indepedent sites within the GMC. Along with other examples, such as the G305 and Carina star-forming regions, we suggest that W51 is a Galactic analogue of the ubiquitous star cluster complexes seen in external galaxies such as M51 and NGC2403.Peer reviewe

Microcassiope minor (Dana, 1852): a description of the first stage zoea (Crustacea: Decapoda: Brachyura: Xanthidae)

The first stage zoea of Microcassiope minor (Dana, 1852) is described, illustrated and compared with thirteen other known Xanthinae. This comparison suggests that M. minor is closely related to another Atlantic Ocean xanthid, Nanocassiope melanodactyla (A. Milne-Edwards, 1867). Such an affinity may indicate that the genus Nanocassiope Guinot, 1967) is not a homogenous taxon

Solving Lattice QCD systems of equations using mixed precision solvers on GPUs

Modern graphics hardware is designed for highly parallel numerical tasks and promises significant cost and performance benefits for many scientific applications. One such application is lattice quantum chromodyamics (lattice QCD), where the main computational challenge is to efficiently solve the discretized Dirac equation in the presence of an SU(3) gauge field. Using NVIDIA's CUDA platform we have implemented a Wilson-Dirac sparse matrix-vector product that performs at up to 40 Gflops, 135 Gflops and 212 Gflops for double, single and half precision respectively on NVIDIA's GeForce GTX 280 GPU. We have developed a new mixed precision approach for Krylov solvers using reliable updates which allows for full double precision accuracy while using only single or half precision arithmetic for the bulk of the computation. The resulting BiCGstab and CG solvers run in excess of 100 Gflops and, in terms of iterations until convergence, perform better than the usual defect-correction approach for mixed precision.Comment: 30 pages, 7 figure

A new dromaeosaurid (Dinosauria: Theropoda) from the Upper Cretaceous Wulansuhai Formation of Inner Mongolia, China

We describe a new dromaeosaurid theropod from the Upper Cretaceous Wulansuhai Formation of Bayan Mandahu, Inner Mongolia. The new taxon, Linheraptor exquisitus gen. et sp. nov., is based on an exceptionally well-preserved, nearly complete skeleton. This specimen represents the fifth dromaeosaurid taxon recovered from the Upper Cretaceous Djadokhta Formation and its laterally equivalent strata, which include the Wulansuhai Formation, and adds to the known diversity of Late Cretaceous dromaeosaurids. Linheraptor exquisitus closely resembles the recently reported Tsaagan mangas. Uniquely among dromaeosaurids, the two taxa share a large, anteriorly located maxillary fenestra and a contact between the jugal and the squamosal that excludes the postorbital from the infratemporal fenestra. These features suggest a sister-taxon relationship between L. exquisitus and T. mangas, which indicates the presence of a unique dromaeosaurid lineage in the Late Cretaceous of Asia. A number of cranial and dental features seen in L. exquisitus and T. mangas, and particularly some postcranial features of L. exquisitus, suggest that these two taxa are probably intermediate in systematic position between known basal and derived dromaeosaurids. The discovery of Linheraptor exquisitus is thus important for understanding the evolution of some salient features seen in the derived dromaeosaurids

Experimental control of a cupola furnace

In this paper the authors present some final results from a research project focused on introducing automatic control to the operation of cupola iron furnaces. The main aim of this research is to improve the operational efficiency and performance of the cupola furnace, an important foundry process used to melt iron. Previous papers have described the development of appropriate control system architectures for the cupola. In this paper experimental data is used to calibrate the model, which is taken as a first-order multivariable system with time delay. Then relative gain analysis is used to select loop pairings to be used in a multiloop controller. The resulting controller pairs melt rate with blast volume, iron temperature with oxygen addition, and carbon composition with metal-to-coke ratio. Special (nonlinear) filters are used to compute melt rate from actual scale readings of the amount of iron produced and to smooth the temperature measurement. The temperature and melt rate loops use single-loop PI control. The composition loop uses a Smith predictor to discount the deadtime associated with mass transport through the furnace. Experiments conducted at the Department of Energy Albany Research Center`s experimental research cupola validate the conceptual controller design and provide proof-of-concept of the idea of controlling a foundry cupola